QUENCHED LIMITS AND FLUCTUATIONS OF THE EMPIRICAL 
MEASURE FOR PLANE ROTATORS IN RANDOM MEDIA 

ERIC LUgON 

Abstract. The Kuramoto model has been introduced to describe synchronization phe- 
nomena observed in groups of cells, individuals, circuits, etc. The model consists of A'^ 
interacting oscillators on the one dimensional sphere S^, driven by independent Brownian 
Motions with constant drift chosen at random. This quenched disorder is chosen inde- 
pendently for each oscillator according to the same law /i. The behaviour of the system 
("^ , for large A^ can be understood via its empirical measure: we prove here the convergence 

Cn ' as TV — >■ oo of the quenched empirical measure to the unique solution of coupled McKean- 

r^ ' Vlasov equations, under weak assumptions on the disorder n and general hypotheses on 

C^ , the interaction. The main purpose of this work is to address the issue of quenched fluctu- 

ations around this limit, motivated by the dynamical properties of the disordered system 
for large but fixed A^: hence, the main result of this paper is a quenched Central Limit 
CO ■ Theorem for the empirical measure. Whereas we observe a self-averaging for the law 

of large numbers, this no longer holds for the corresponding central limit theorem: the 
J , trajectories of the fluctuations process are sample-dependent. 

1. Introduction 

In this work, we study the fluctuations in the Kuramoto model, which is a particular 
Cn ' case of interacting diffusions with a mean fleld Hamiltonian that depends on a random 

disorder. The Kuramoto model was first introduced in the 70's by Yoshiki Kuramoto 
^T) [ ([17j. see also [T] and references therein) to describe the phenomenon of synchronization 

C^ ' in biological or physical systems. More precisely, the Kuramoto model is a particular 

case of a system of N oscillators (considered as elements of the one-dimensional sphere 
P: ; S^ := R/27rZ) solutions to the following SDE: 



O 



X 



1 ^ 
(1) dxl'^ = —^bix'^^,x^^^,ujj)dt + cix'^^,uJi)dt+dBi, t£[0,T], i = l...N, 



N 



where T > is a fixed (but arbitrary) time, b and c are smooth periodic functions. The 
5^ I Kuramoto case corresponds to a sine interaction (b{x,y,uj) = Ksm{y—x) and c{x,uj) = co). 

This case which has the particularity of being rotationally invariant (namely, if (j;*'^)j is 
a solution of the Kuramoto model, (x*' + c)j, c a constant, is also a solution), will be 
referred to in this work as the sine-model. 

The parameter i^ > is the coupling strength and {ujj) is a sequence of randomly 
chosen reals (being i.i.d. realizations of a law ^u). The sequence {ujj)j is called disorder 
and represents the fact that the behaviour of each rotator x^'^ depends on its own local 
frequency coj. 

Due to the mean field character of ([T]), the behaviour of the system can be understood 
via its empirical measure u , process with values in A^i(S"'^ x R), that is the set of 
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probability measures on oscillators and disorder: 



1 ^ 
V(.;) G R^, Vt G [0,r], .f'^^) := ^ E'^C.--,.,)' 

i=i 

where (w) = (wj)j>i is a fixed sequence of disorder in R'^. 

The purpose of this paper is to address the issue of both convergence and fluctuations of 
the empirical measure, as A^ — )■ cxd ; thus the main theorem of this paper (Theorem 12. lOp 
concerns a Central Limit Theorem in a quenched set-up (namely the quenched fluctuations 
of V around its limit). 

Some heuristic results have been obtained in the physical literature ([1] and references 
therein) concerning the convergence of the empirical measure, as A^ — )■ cx), to a time- 
dependent measure {Pt{dx^ da;))(g[o.j'], whose density w.r.t. Lebesgue measure at time i, 
qt{x,uj) is the solution of a deterministic non-linear McKean-Vlasov equation (see Eq.([5])). 
It is well understood ( [Ij , [H] ) that crucial features of this equation are captured in the 
sine-model by order parameters rt and -i/'t defined by: 

rje*'^' = / e^^qt{x,uj)dxn{du}). 

The quantity rt captures in fact the degree of synchronization of a solution (the profile 
qt = -^ for example corresponds to r = and represents a total lack of synchronization) 
and ipt identifies the center of synchronization: this is true and rather intuitive for unimodal 
profiles. Moreover ([22|, see also pUj, p. 118) it turns out that if // is symmetric, all the 
stationary solutions can be parameterized (up to rotation) by the stationary version r of 
rt which must satisfy a fixed point relation r = ^/^^^(r), with ^ k,^{') an explicit function 
such that ^;^^^(0) = 0. For K small, r = is the only solution of such an equation and the 
system is not synchronized, but for K large, non-trivial solutions appear (synchronization). 
In the easiest instances, such a non-trivial solution is unique (in the sense that r = ^ k,^{t~) 
admits a unique non-zero solution but of course one obtains an infinite number of solutions 
by rotation invariance that is ijj can be chosen arbitrarily). 

In [9], Dai Pra and den Hollander have rigorously shown the convergence of the aver- 
aged empirical measure L G A^i(C([0, T], S"*^) x R) (probability measure on the whole 
trajectories and the disorder): 



L^ 



1 ^ 



This convergence of the law of L under the joint law of the oscillators and the disorder is 
shown via an averaged large deviations principle in the case where b{x, y,Lo) = K- f{y— x) 
and c{x,Lo) = g{x,u}) for / and g smooth and bounded functions. As a corollary, it 
is deduced in [9] the convergence of Ljv and of u , via a contraction principle, in the 
averaged set-up. In the case of unbounded disorder, the same proof can be generalized 
(thesis in preparation) under the following assumption: 

(H^) Vi>0, /" e*l'^l/x(da;) <oo. 

One aim of this paper is to obtain the limit of v^'^'^' in the quenched model, namely for 
a, fixed realization of the disorder [uj). This result can be deduced from the large deviations 
estimates in [9j, via a Borel-Cantelli argument, but our result is more direct and works 
under the much weaker assumption on ^u, J-^ |a;|;u(da;) < oo. 
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Figure 1. We plot here the evolution of the marginal on S^ of v^'^^i 
for A^ = 600 oscillators in the sine-model (^ = \{5-\ + (5i), K = 6). 
The oscillators are initially chosen independently and uniformly on [0, 2ii\ 
independently of the disorder. First the dynamics leads to synchronization 
of the oscillators (t = 6) to a profile which is close to a non-trivial stationary 
solution of McKean-Vlasov equation. We then observe that the center 
■04 of this density moves to the right with an approximately constant 
speed; what is more, this speed of fluctuation turns out to be sample- 
dependent (see Fig. [2]). 

A crucial aspect of the quenched convergence result, which is a law of large numbers, is 
that it shows the self- averaging character of this limit: every typical disorder configuration 
leads as A^ — )• oo to the same evolution equation. 

However, it seems quite clear even at a superficial level that if we consider the central 
limit theorem associated to this convergence, self-averaging does not hold since the fluctu- 
ations of the disorder compete with the dynamical fluctuations. This leads for example to 
a remarkable phenomenon (pointed out e.g. in [2] on the basis of numerical simulations): 
even if the distribution /i is symmetric, the fluctuations of a fixed chosen sample of the 
disorder makes it not symmetric and thus the center of the synchronization of the system 
slowly (i.e. with a speed of order 1/^/N) rotates in one direction and with a speed that 
depends on the sample of the disorder (Fig. [T]and[2]). This non-self averaging phenomenon 
can be tackled in the sine- model by computing the finite-size order parameters (Fig. [2]): 



1 



N 

N 

As a step toward understanding this non self-averaging phenomenon, the second and 
main goal of this paper is to establish a fluctuations result around the McKean-Vlasov 
limit in the quenched set-up (see Theorem I2.10|) . A central limit theorem for the averaged 
model is addressed in [S] , applying techniques introduced by Bolthausen [B] . A fluctuations 
theorem may also be found in [S] for a model of social interaction in an averaged set-up. 
Here we prove convergence in law of the quenched fluctuations process 



r?^.H:=Viv(i.^'(-)-P) 



ERIC LUgON 




Figure 2. Trajectories of the center of synchronization ijj"'^^) in the sine- 
model for different reaUzations of the disorder (/i = i((^_o.5 + ^0.5), /f = 4, 
A'^ = 400). We observe here the non self-averaging phenomenon: direction 

and speed of the center depend on the choice of the initial A'^-sample of the 
disorder. Moreover, these simulations are compatible with speeds of order 

1/VN. 



seen as a continuous process in the Schwartz space S' of tempered distributions on S^ x R 
to the solution of an Ornstein-Uhlenbeck process. The quenched convergence is here un- 
derstood as a weak convergence in law w.r.t. the disorder and is more technically involved 
than the convergence in the averaged system. The main techniques we exploit have been 
introduced by Fernandez and Meleard |l12j and Hitsuda and Mitoma ^15j, who studied 
similar fluctuations in the case without disorder. In [7], A Large Deviation Principle is 
also proved. We refer to Section H] for detailed definitions. 

While numerical computations of the trajectories of the limit process of fluctuations 
clearly show a non self-averaging phenomenon, the dynamical properties of the fluctuations 
process that we find are not completely understood so far. Progress in this direction 
requires a good understanding of the spectral properties of the linearized operator of 
McKean-Vlasov equation around its non-trivial stationary solution; the stability of the 
non-synchronized solution q = ^ has been treated by Strogatz and Mirollo in [25 . In the 
particular case without disorder, spectral properties of the evolution operator linearized 
around the non-trivial stationary solution are obtained in [3J, but the case with a general 
distribution /x needs further investigations. 

This work is organized as follows: Section [2] introduces the model and the main results. 
Section [3] focuses on the quenched convergence of i^jy- In Section HI the quenched Central 
Limit Theorem is proved. The last section applies the fluctuations result to the behaviour 
of the order parameters in the sine-model. 



2. Notations and main results 

2.1. Notations. 

• if X is a metric space, Bx will be its Borel cr-field, 



QUENCHED LIMITS AND FLUCTUATIONS FOR ROTATORS IN RANDOM MEDIA 



• Ch{X) (resp. C^{X), p = 1, . . . , cx)), the set of bounded continuous functions (resp. 
bounded continuous with bounded continuous derivatives up to order p) on X, {X 
will be often S^ x R), 

• Cc{X) (resp. Cc{X), p = 1, . . . ,00), the set of continuous functions with compact 
support (resp. continuous with compact support with continuous derivatives up 
to order p) on X, 

• D([0,T],X), the set of right-continuous with left limits functions with values on 
X, endowed with the Skorokhod topology, 

• Mi{Y), the set of probability measures on Y (Y topological space, with a regular 
a-fieldi3), 

• Mf{Y), the set of finite measures on Y, 

• {Aii{Y),w): Aii(Y) endowed with the topology of weak convergence, namely the 
coarsest topology on M.i{Y) such that the evaluations v h^ J f du are continuous, 
where / are bounded continuous, 

• {Aii{Y),v): M.i{Y) endowed with the topology of vague convergence, namely the 
coarsest topology on M.i{Y) such that the evaluations i' ^ j f du are continuous, 
where / are continuous with compact support. 

We will use C as a constant which may change from a line to another. 

2.2. The model. We consider the solutions of the following system of SDEs: 
for z = 1, . . . , iV, for T > 0, for ah t < T, 

-I ^ ft ft 

(3) xj'^ = ^^ + -Y,l b(4'^, ^i'^, ^j) ds + j^ c(xi'^, u;,) ds + Bl 



3-- 

where the initial conditions i^ are independent and identically distributed with law A, and 
independent of the Brownian motion (5) = (i?')j>i, and where h (resp. c) is a smooth 
function, 27r-periodic w.r.t. the two first (resp. first) variables. The disorder {uS) = (ciJj)j>i 
is a realization of i.i.d. random variables with law ^i. 
Remark 2.1. 

The assumption that the random variables (wj) are independent will not always be neces- 
sary and will be weakened when possible. 

Instead of considering x*' as elements of R, we will consider their projection on S^. 
For simplicity, we will keep the same notation x*' for this projectiorU- 

We introduce the empirical measure v (on the trajectories and disorder): 
Definition 2.2. 

For all t <T, for a fixed trajectory {x^, . . . , x^) G C([0, T], (S^)^) and a fixed sequence of 
disorder {co), we define an element of A4i{S^ x R) by: 

Finally, we introduce the fiuctuations process rj^'^'^' of u^'^'^' around its limit P (see 
Th.ESD: 
Definition 2.3. 
For all t <T, for fixed (uj) G R , we define: 

Throughout this article, we will denote as P the law of the sequence of Brownian Motions 
and as P the law of the sequence of the disorder. The corresponding expectations will be 
denoted as E and E respectively. 



See Remarks 12. 71 and 12. 121 for possible generalizations to the non-compact case. 



ERIC LUgON 



2.3. Main results. 

2.3.1. Quenched convergence of the empirical measure. In [5], Dai Pra and den Hollander 
are interested in the averaged model, i.e. in the convergence in law of the empirical 
measure under the joint law of both oscillators and disorder. The model studied here, 
which is more interesting as far as the biological applications are concerned is quenched: 
for a fixed realization of the disorder (to), do we have the convergence of the empirical 
measure? Moreover the convergence is shown under weaker assumptions on the moments 
of the disorder. 

We consider here the general case where b(x, y, uj) is bounded, Lipschitz-continuous, and 
27r-periodic w.r.t. the two first variables, c is assumed to be Lipschitz-continuous w.r.t. 
its first variable, but not necessarily bounded (see the sine-model, where c{x,ij) = uj). 
We also suppose that the function oj i— ;• S{ijj) := supj,gsi |c(x,6t;)| is continuous (this is in 
particular true if c is uniformly continuous w.r.t. to both variables {x,uj), and obvious in 
the sine-model where S{uj) = \uj\). The Lipschitz bounds for b and c are supposed to be 
uniform in oj. 

The disorder (w), is assumed to be a sequence of identically distributed random variables 
(but not necessarily independent), such that the law of each w, is //. We suppose that the 
sequence (a;) satisfies the following property: for P-almost every sequence (w), 

1 ^ f 

(H^) T7y~! sup |c(x,Wj)| ->Af-_>oo / sup \c{x,uj)\fi{duj) < oo. 

We make the following hypothesis on the initial empirical measure: 

(Hq) Uq — > vq, in law, in (A^i(S x R),t(;). 

Remark 2.4. (1) The required hypotheses about the disorder and the initial conditions 
are weaker than for the large deviation principle: 

• the (identically distributed) variables (wj) need not be independent: we simply 
need a convergence (similar to a law of large numbers) only concerning the 
function S, 



Condition (H/^ ) is weaker than {H^) on page EJ for the sine- model, {H^) 



reduces to f \uj\fi{duj) < oo, 

the initial values need not be independent, we only assume a convergence of 

the empiric al m easure. 



(2) The hypothesis (Hj2 ) is verified, for example, in the case of i.i.d. random variables, 
or in the case of an ergodic stationary Markov process. 

(3) Under ( |ffoD , the second marginal of vq is fi. 



In Section [3l we show the following: 



Theorem 2.5. 



Under the hypothesis l\Ho\ i and {HJj), forF-almost every sequence (uji), the random vari- 



able i/^'^^i converges in law to P, in the space D([0,T], (A^i(S^ x R),i(;)), where P is the 
only solution of the following weak equation (for every f continuous bounded on S"*^ x R, 
twice differentiable, with bounded derivatives): 

(4) {Pt ,f) = {uo,f) + l[ ds (P, , f") + f ds{Ps, nb[;Ps] + C)), 

■^ Jo Jo 



where 



b[x,m] = / b{x,y,Tr)m{dy, dvr). 
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Moreover, with the same hypotheses, the law of v^ under the joint law of the oscillators 
and disorder (averaged model) converges weakly to P as well. 

Remark 2.6. 

An easy calculation shows that P can be considered as a weak solution to the family of 

coupled McKean-Vlasov equations (see [9]): 

(1) P can be written as P(dx, do;) = //(a;)P'^(dx), 

(2) if we define q^ through Pt{dx, dw) = /x(w)g^(dx), q^ is the unique weak solution 
of the McKean-Vlasov equation: 

(5) ^^qt = C-q^. q^ = X. 



where, C^ is the following differential operator: 

b{x, y, ■^)qt ( dy)/x( dvr) + c(x, uj) ] qf 



2 



1 d' 



We insist on the fact that Eq. ^ is indeed a (possibly) infinite system of coupled non- 
linear PDEs. To fix ideas, one may consider the simple case where ^u = 2('^-i + "^i)- Then 
([5]) reduces to two equations (one for -|-1, the other for —1) which are coupled via the 
averaged measure \{qf + Qt )• ^^^ fo'^ more general situations (/i = AA(0, 1) say) this 
would consist of an infinite number of coupled equations. 

Remark 2.1 (Generalization to the non compact case). 

The assumption that the state variables are in S^, although motivated by the Kuramoto 
model, is not absolutely essential: Theorem 12.51 still holds in the non-compact case (e.g. 
when S^ is replaced by R'^), under the additional assumptions of boundedness of x i— )• 
|c(x,a;)| and the first finite moment of the initial condition: J |x|A(dx) < oo. 

We know turn to the statement of the main result of the paper: Theorem 12.101 

2.3.2. Quenched fluctuations of the empirical measure. Theorem [23] says that for P-almost 
every realization (w) of the disorder, we have the convergence of u^'^'^' towards P, which 
is a law of large numbers. We are now interested in the corresponding Central Limit 
Theorem associated to this convergence, namely, for a, fixed realization of the disorder (w), 
in the asymptotic behaviour, as iV -^ oo of the fluctuations field rj^'^^i taking values in 
the set of signed measures: 

Vt G [0,r], r?f •(") := Vn fi/f'^") - P, 



In the case with no disorder, such fluctuations have already been studied by numer- 
ous authors (eg. Sznitman [25], Fernandez- Meleard [12], Hitsuda-Mitoma [15j). More 
particularly, Fernandez and Meleard show the convergence of the fluctuations field in an 
appropriate Sobolev space to an Ornstein-Uhlenbeck process. 

Here, we are interested in the quenched fluctuations, in the sense that the fluctuations 
are studied for fixed realizations of the disorder. We will prove a weak convergence of the 
law of the process r/^'^'^), in law w.r.t. the disorder. 

In addition to the hypothesis made in ^2.3. H we make the following assumptions about 
h and c (where Vp is the set of all differential operators of the form d^kd^i with k + l < p): 



iHb,c) 



beC^iS^ xK), cgC~(S1xR), 



DeVeJR 1 + kl 



3a > 0, sup / -, , i_|2a' dvr < oo, 
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Furthermore, we make the fohowing assumption about the law of the disorder (a is defined 



in (Fb,j): 



(//^) the {loj) are i.i.d. and / |a;| "^(do;) < qg. 

Jr. 

Remark 2.8. 

The regularity hypothesis about b and c can be weakened (namely b G C^(S^ x R) and 

c G ^"(S^ X R) for sufficiently large n and m) but we have kept m = n = oo ioi the sake 

of clarity. 

Remark 2.9. 



In the case of the sine-model, Hypothesis (H^c) is satisfied with a = 2 for example. 



In order to state the fluctuations theorem, we need some further notations: for all s <T, 
let Cs be the second order differential operator defined by 

^s{^){y,-^) ■■= -(p"{y,TT) + ip'{y,'7T){b[y,Ps] + c(y,7r)) + {Pg , ip' {■ , ■)b{- , y , it)) . 
Let W the Gaussian process with covariance: 

(7) B{Wt{vi)Ws{v2)) = {Pu,vW2)du. 

Jo 

For all (/?!, ip2 bounded and continuous on S^ x R, let 

(8) Ti{ipi,ip2)= Covxifii- ,u;),(p2i- ,uj))fi{du;), 

Jn 

((/?i-/ (pi{- ,uj)dX] [(f2- (f2{- ,^)<iX] X{dx)fi{duj), 
SixR V ^si / V ^si / 



and 

(9) r2{(pi,(f2) = Cov^ I / (^idA, / ip2dX] 



V?! dA — / cpidXdfi ] I ip2dX— c^2dAd/i) dfi. 

R \Js^ JsixR / VJsi JsixR / 

For fixed (cj), we may consider ^{^{uj), the law of the process 7/^'('^); %]^[uj) belongs 
to 7Wi(C([0,r],5')), where S' is the usual Schwartz space of tempered distributions on 
S^ X R. We are here interested in the law of the random variable (w) i— )■ T-l]\f{uj) which is 
hence an element of Mi{Mi{C{[0,T],S'))). 

The main theorem (which is proved in Section [3|) is the following: 

Theor em 2 .10 (Fluctuations in the quenched model). 



Under (H^), (i^b^d, the sequence (uj) i— )• ^^{ijj) converges in law to the random variable 
uj I— 7- T-L{uj), where l-L{oj) is the law of the solution to the Ornstein-Uhlenbeck process r]^ 
solution in S' of the following equation: 

(10) rjf=Xiu;)+ f £:<ds + T^t, 

Jo 

where, £* is the formal adjoint operator of Cs and for all fixed u, X{ui) is a non-centered 

Gaussian process with covariance Fi and with mean value C{uj). As a random variable in 

CO, CO >-^ C{u}) is a Gaussian process with covariance T2. Moreover, W is independent on 

the initial value X . 

Remark 2.11. 

In the evolution (|10p . the linear operator £* is deterministic ; the only dependence in oj 
lies in the initial condition X{(jj), through its non trivial means C{uj). However, numerical 
simulations of trajectories of r/"^ (see Fig. [3]) clearly show a non self-averaging phenomenon, 
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analogous to the one observed in Fig [2j Vj^ not only depends on uj through its initial 
condition X{uj), but also for all positive time t > 0. 

Understanding how the deterministic operator >C* propagates the initial dependence in 
uj on the whole trajectory is an intriguing question which requires further investigations 
(work in progress). In that sense, one would like to have a precise understanding of the 
spectral properties of £*, which appears to be deeply linked to the differential operator in 
McKean-Vlasov equation ([6]) linearized around its non-trivial stationary solution. 

Remark 2.12 (Generalization to the non-compact case). 

As in Remark 12.71 it is possible to extend Theorem 12.101 to the (analogous but more 
technical) case where S^ is replaced by R . To this purpose, one has to introduce an 
additional weight (1 + |rE|°)~^ in the definition of the Sobolev norms in Section [Hand to 
suppose appropriate hypothesis concerning the first moments of the initial condition A 
(J |x|'^A(dx) < oo for a sufficiently large /3). 

2.3.3. Fluctuations of the order parameters in the Kuramoto model. For given A^ > 1, 
t G [0,T] and disorder (w) G R^, let rf'^""^ > and cf'^""^ G S^ such that 



N 



Af,H>Af,(a;) _ I sr^ j^ 



Proposition 2.13 (Convergence and fluctuations for r^ ). 

We have the following: 

(1) Convergence of r^ : For¥-almost every realization of the disorder, r '^''^> con- 
verges in law in C([0,T],R), to r defined by 

1 
t G [0,r] ^ rt := ((Pi , cos(-))' + {Pt , sin(-))'^ ' 

(2) // ro > then 
(Hr) VtG[0,r], rt>0. 

(3) Fluctuations of r^ around its limit: Let 

t \-^TZ^ := vN {r^ — rt 

be the fluctuations process. For fixed disorder (w), let 9^^'('^) £ A//i(C([0,T],R)) 
be the law of TZ '^'^' . Then, under (\Hr\ i, the random variable [uj) i— )• y{"'^^> con- 
verges in law to the random variable uj i— )• 9^'^, where yi'^ is the law of IZ^ := 
i {{P, cos(-)) • (r/- , cos(-)) + {P, sin(-)) • « , sin(.))). 

Remark 2.14- 

In simpler terms, this double convergence in law corresponds for example to the conver- 
gence in law of the corresponding characteristic functions (since the tightness is a direct 
consequence of the tightness of the process r/); i.e. for ti,...,tp £ [0,T] (p > 1) the 

characteristic function of (7^^ ^ , . . . ,7^j ) for fixed (a;) converges in law, as a random 
variable in (oj), to the random characteristic function of {TZf^, . . . ,TZf ). 

Proposition 2.15 (Convergence and fluctuations for ^^'(^)). 

We have the following: 

(1) Convergence of C, '^'^' : Under ^Hj.\ , for¥-almost every realization of the disorder 
(uj), ^"'^'^^ converges in law to Q : t £ [0,T] i— t- C,t '■= — - — , 



10 
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>(2r; 




15 
t 



Figure 3. We plot here the evolution of the imaginary part of the process 
Z^, for different realizations of uj; the trajectories are sample-dependent, 
as in Fig [21 

(2) Fluctuations of C^'('^) around its limit: Let 



t^ Z. 



N,{uj) 



iv(cf'(-)-G) 



be the fluctuations process. For fixed disorder (uj), let 3^'^'^'' € 7Wi(C([0, r],R)) 
be the law of Z^'^'^K Then, under ( |ff^[ ), the random variable (a;) i— t- 3^''"^^ con- 
verges in law to the random variable uj ^^ '^'^ , where 3"^ is the law of Z"^ := 
;^(r(r?-,cos(•)) + <P,e">7^-). 

In the sine- model, we have Cj^-'^^i = e^^ ' " where ijj^'^'^' is defined in ([2]) and is plotted 
in Fig. [2j Some trajectories of the process Z^ are plotted in Fig. [3l 
This fluctuations result is proved in Section [5l 

3. Proof of the quenched convergence result 
In this section we prove Theorem 12.51 Reformulating ([3]) in terms of u '^^\ we have: 



(11) vi = i...iv, vtG [o,r], 



i,N 



e + 1 b[x 

'0 



i,N .,N 



, i/g ] ds + / c[x 





i,N 



,u}i)ds + Bi, 



where we recall that 6[x,rn-] := J b{x,y,TT)m{dy, dvr). 

The idea of the proof of Theorem 12.51 is the following: we show the tightness of the se- 
quence (u '^^') firstly in D([0, T], {AiF,v)) (recall Notations in §2. ID . which is quite simple 
since Cc(S^ x R) is separable and by a n arg ument of boundedness of the second marginal 
of any accumulation point, thanks to (Hjj), we show the tightness in D([0,r], {M-f,w)). 



The proof is complete when we prove the uniqueness of any accumulation point. 

3.1. Proof of the tightness result. We want to show successively: 

(1) Tightness of £(z^^'(")) in B{[0,T],{Mf,v)), 

(2) Equation verified by any accumulation point, 

(3) Characterization of the marginals of any limit, 

(4) Convergence in D([0,r], {Mf,w)). 
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3.1.1. Equation verified by i/^'('^). For / G C^(S^ x R), we denote by /', /" the first and 
second derivative of / with respect to the first variable. Moreover, if m £ A^i(S^ x R), 
then {m , f) stands for J^^^-^ f{x,TT)m{dx, dvr). 

Applying Ito's formula to (fTTI) . we get, for all / G C^(S^ x R), 



Af,(a;) 



/ 



Af,(a;) 



/ + 



\ ['<""■'"'■ f) 



+ 



ds (z.f 'H , /' • {b[; Z.f 'H] + C)) + M^jit), 





rt fu N,{uj) 



where Mj^j{t) := jj X^j^^ /o f'i^j '"^i) d-Bj(s) is a martingale (/' bounded). 

3.1.2. Tightness of £{iy^'^'^'^) in B{[0,T],{Mf,v)). Cc(S^ x R) is separable: let {fk)k>i 
(elements ofC°°(S^ ^R-)) a dense sequence in Cc(S^ x R), and let /o = 1. We define 
n := D([0,r], {Mi,v)) and the applications Uj, f £ Cc{S^ x R) by: 



Uf : Q 



D([0,r],R) 



m I— 7- {m , f) . 

Let {Pn)n a sequence of probabilities on il and (HfPn) = Pn ° 117 G D([0,T],R). We 
recall the following result: 

Lemma 3.1. 

If for all k > 0, the sequence (njj.P„)„ is tight in A^i(D([0, T], R)), then the sequence 
{Pn)n is tight m Mi{T){[id,Tl{Mi,v))). 

Hence, it suffices to have a criterion for tightness in D([0, T], R). Let X" be a sequence 
of processes in D([0, T], R) and J^" a sequence of filtrations such that X" is J^^-adapted. 
Let (/)" = {stopping times for J-""}. We have (cf. Billingsley ^): 

Lemma 3.2 (Aldous' criterion). 

// the following holds, 

(1) C (supt<r |Xf 1)^ is tight, 

(2) For all e > and r] > 0, there exists 5 > such that 

limsup sup P {\X^ -X^,\>rj)< e, 

n 5,5'6(/>";S<S'<(5+<5)AT 
then C{X"') is tight. 

Proposition 3.3. 

The sequence £(z^^'('^)) is tight in D([0,r], {Mf,v)). 

Proof. For all e > 0, for all /c > 1 (the case k = is straightforward). 



P sup 

\t<T 



^t ' Jk) > ~ ] ^ ^W Jk 



E 



sup 

t<T 



.f ■(-) , 1 



, (Markov Inequality). 



The tightness follows. 

For all A; > 1, we have the following decomposition: 



N,{ui) 



fk 



'^0 ^ fk) + A Uk) + Mt '{fk), 



,7V, (a;) 



NMl 



where A^ (fk) is a process of bounded variations, and M^ (fk) is a square-integrable 
martingale. Then it suffices to verify Lemma 13.21 ^ for A and M separately. For all 
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£ > and 77 > 0, for all stopping times S, S' G </>^; S < S' < {S + 6) AT, we have: 
a„:=p(|<;<"'(A)-Af"'(A)|>,). 



<iE 

V 



s' 



ds 



^f'(^^/^(&h^f'(")]+c) 






S' 



ds 



"^s ^ Jk 



< — E \S' — S] < e, for 6 sufficiently small. 

7] '- -' 

(we use here that fk are of compact support for A; > 1; in particular the function (x, vr) 1— >• 
f'f^{x,iT)c{x,iT) is bounded). Furthermore, 



1 



<-E 



< 



{Nr,y 



rE 



~1 Js 



12 fi 



x\,ijJi) ds 



> rj 

2" 



<-^d. n 



Nt]^ 



At this point, £(zy^'H) is tight in D([0,r], {Mf,v)). 



3.1.3. Equation satisfied by any accumulation point in D([0,T], {Mf,v)). Using hypoth- 
esis dffoD , it is easy to show that the following equation is satisfied for every accumulation 
point u, for every / G C^i^^ x R) (we use here that S^ is compact): 

(12) {ut ,f) = {uo,f)+f ds {us , /' • (6[-, vs\ + c)) + \ f ds {v, , /">. 

JO ^ Jo 



For any accumulation point u, the following lemma gives a uniform bound for the second 
marginal of i^: 

Lemma 3.4. 

Let Q be an accumulation point of C{v^'^^'>)n in D([0,T], {M.i,v)) and let be v ^ Q. For 
all t € [0, r], we define by (ft, 2) ihe second marginal of Vf' 



yAeBiK), {vt,2){A)= [ vt{dx,dn). 



Then, for all t G [0, T], 



sup |c(x,7r)| (i/i,2)(d7r) < / sup |c(x,7r)| /i(d7r). 



Proof. Let (/> be a C^ positive function such that cp = 1 on [—1,1], = on [—2,2] and 
U\\oo<'^- Let, 



VA; > 1, (j)k ■■= IT ^ (/) (-] . 



Then (pk G C^(R-) and </>fe(7r) 



^fc— >oo 



1, for all vr. We have also for all vr G R, |(/>fc(7r)| < 



^\L^m^)\<u'\\oo^m^)\<u" 
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We have successively, denoting ^'(vr) := sup^ggi |c(x,7r) 



S{iT)i/t{dx, dn) = / limini (f)k{T^)S{-iT)i't{dx, dir), 

SixR JSixR k^'^ 

< liminf / (f)k{T^)S{-iT)i't{dx, dir), (Fatou's lemma), 

fc^oo JgixR 

(13) = liminf lim / (j)k{'K)S{'K)v^ (dx, dvr), 

< lim / S{'k)vi "^ (dx, dvr), (since ||</'||oo < !)• 

Af^ooJgixR, 

The equality (fT3]l is true since (x, vr) i— )■ 0fc(7r)5(7r) is of compact support in S^ x R (recall 
that S is supposed to be continuous by hypothesis). 



But, by definition of v^ , and using the hypothesis (Hji) concerning //, we have. 



(14) lim / 5(7r)z.f'('^^(dx, d7r)= / S(^)/u(d^). 

The result follows. D 

Remark 3.5. 

(fH|l is only true for P- almost every sequence (uj). We assume that the sequence (w) given 

at the beginning satisfies this property. 

3.1.4. Tightness in D([0,T], {Mf,w)). We have the following lemma (cf. (21]): 

Lemma 3.6. 

Let (Xn) be a sequence of processes in DdOjT], {A4p,w)) and X a process belonging to 
C{[Q,T],{Mf,w)). Then, 

X^X^i^^^^ mD([0,r],(A^^,t;)), 

1 {Xn, 1)^{X, 1) mD([0,r],R). 

So, it suffices to show, for any accumulation point z^: 

(1) {v , 1) = 1: Eq. ([12]) is true for all / G C^iS^ x R), so in particular for fk{x,iT) := 
(j)k{TT). Using the boundedness shown in lemma 13. 4| we can apply dominated 
convergence theorem in Eq. (I12p . We then have {vt , 1) = 1, for all t £ [0,T]. The 
fact that Eq. (I12p is verified for all / G C^(S^ x R) can be shown in the same way. 

(2) Continuity of the limit: For ah < s < t < T, for all / G C^(S^ x R), 

1 '■' 



\{i^t,f)-{i^s,f)\<K J \{i,^,f'.b[-,i^u])\du+-J \{iy^,f")\du 

+ / {^u ) /' ■ c) du < C X |t — s| , for some constant C. 
J s 

Noticing that we used again Lemma [3.4l to bound the last term, we have the result. 

Remark 3. 7. 

It is then easy to see that the second marginal (on the disorder) of any accumulation point 

is /i. 

At this point C{v^'^'^>) is tight in D([0, T], (7W_p, w)). It remains to show the uniqueness 
of any accumulation point: it shows firstly that the sequence effectively converges and that 
the limit does not depend on the given sequence (w). 
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3.2. Uniqueness of the limit. 
Proposition 3.8. 

There exists a unique element P o/D([0,T], A^i7'(S^ x R)) which satisfies equation (J12p . 
Pq e Ml and Po,2 = A*- 

Remark 3.9. 

Since this proof is an adaptation of Oelschlager |20j Lemma 10, p. 474, we only sketch the 

proof: 

We can rewrite Eq. ()12p in a more compact way: 



(15) {Pt ,f) = {Po,f)+[ (P, , L{Ps){f)) ds, V/ G C^S' X R), < t < T, 

Jo 



where. 



L{P){f){x,uj) := ^f"{x,oj) + h{x,io,P)f'{x,u), Vx G S\V^ G R, 



and, 

h{x, vr, P) = b[x, P] + c(x, vr). 
Let t I— )• Pt be any solution of Eq. (|15|) . One can then introduce the following SDE 
(where ^ G R and ^ G S^ is its projection on S^): 

(^f.^ j d^t = h{it,uJt,Pt)dt+ dWt, eo = fo, (lo,^o)~i^o, 

^'""^ \ dcjt = 0. 

Eq. (fT6]l has a unique (strong) solution (Ct)'^t)tG[o,T] = {S,t,^o)te[o,T]- The proof of unique- 
ness in Eq. (|12p consists in two steps: 

(1) Pt = £(6, wt), for ah tG[0,r], 

(2) Uniqueness of ^. 

We refer to Oelschlager [20] for the details. Another proof of uniqueness can be found 
in [9 1 or in [13] via a martingale argument. 

4. Proof of the fluctuations result 



Now we turn to the proof of Theorem 12. 101 To that purpose, we need to introduce some 
distribution spaces: 

4.1. Distribution spaces. Let S := 5(S^ x R) be the usual Schwartz space of rapidly 
decreasing infinitely differentiable functions. Let Dp be the set of all differential operators 
of the form d^kd^i with k + I < p. We know from Gelfand and Vilenkin [_14j p. 82-84, 
that we can introduce on 5 a nuclear Frechet topology by the system of seminorms || • || , 
p = 1,2, . . . , defined by 



/ (l + |vr|2)2p Y, \Dcl,{y,7T)\^dyd7T. 



P 



Let S' be the corresponding dual space of tempered distributions. Although, for the 
sake of simplicity, we will mainly consider rj"'^'^> as a process in C([0, T], 5'), we need some 
more precise estimations to prove tightness and convergence. We need here the following 
norms: 

For every integer j, a G R"^, we consider the space of all real functions (/? defined on 
S^ X R with derivative up to order j such that 

. 1/2 
V^ / \d:j,k^d^k^ip{x,TT)\'^ 



^«^."- ,4,A..^TT#-^-*l <- 
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Let Wq°' be the completion of C^(S^ x R) for this norm; {Wq°', \\ ■ \\j a) is a Hilbert space. 

Let Wq-''"' be its dual space. 

Let C^'"' be the space of functions (p with continuous partial derivatives up to order j 
such that 

hm sup ' ^ "" , , ^' = 0, for all ki + k2 < j, 

kK«=zGSi 1 + KI" 



with norm 

xesiTTsR l + kl" 



V'llci." = X] ^^P ^^P 



ki+k2<j 

We have the following embeddings: 

^m+j> ^ ^j,a^^ > 1, j > 0,a > 0, 

i.e. there exists some constant C such that 

f 17) II UP Wryi r. < C II W II 

V V W r IIGJ'" — ^ W r \\m+j,a 

Moreover, 

^^m+j,a ^ ^,,a+;3^^ ^ 1, J > 0, O > 0, /3 > L 

Thus there exists some constant C such that 

II ^\\j,a+l3 ^ ^11 ^\\m+j,a ■ 

We then have the following dual continuous embedding: 

(18) Tyo--''°+^ ^ W-^"'^'^''', m > 1, a > 0, /3 > 1. 

It is quite clear that 5 M- Wq'" for any j and a, with a continuous injection. 
We now prove some continuity of linear mappings in the corresponding spaces: 
Lemma 4.1. 

For every x,y ^ S^ , uj ^ H, for all a, the linear mappings Wq'°^ — )• R defined by 

Dx,y,u}{v) := ^ix,uj) -ip{y,uj);D^^^ := ip{x,uj); H^^^ = ip'{x,uj), 
are continuous and 

(19) \\D,^y,^\\_^^^<C\x-y\{l + \un, 

(20) ||L>,,^|L3^^<C(1 + |l^|"), 

(21) ||F,,^||_3^„<C(l + |wr). 

Proof. Let c^ be a function of class C°° with compact support on S^ x R, then, 

\ip{x,uj) - ip{y,uj)\ <\x-y\ sup \ip'{u,uj)\ , 

u 

\ip'{u,uj)\' 



< |a; — y| (1 + |a;|") sup 



^ ^iii I) I - 

U,UJ \ 1 + If^l 



< |x-y|(l + |a;|")||(^||i^„, 
<C|x-2/|(l + |^r)||(^||3_^ 



following (I17p with j = 1 and m = 2 > 1. Then, ()19p follows from a density argument. 
(I20p and (I2ip are proved in the same way. D 
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4.2. The non-linear process. The proof of convergence is based on tlie existence of the 
non-hnear process associated to McKean-Vlasov equation. Such existence has been studied 
by numerous authors (eg. Dawson [TO] , Jourdain-Meleard [16] , Mahieu [18] , Shiga- Tanaka 
|23] , Sznitman [25] , [26] ) mostly in order to prove some propagation of chaos properties in 
systems without disorder. We consider the present similar case where disorder is present. 
Let us give some intuition of this process. One can replace the non-linearity in Eq. ([5]) by 
an arbitrary measure m{dx, dcj): 

dtqt = Id^.qt - dx Mx, m] + c{x, u)) q^] . 

In this particular case, it is usual to interpret q^ as the time marginals of the following 
diffusion: 

dx^ = dBt + h[x^ ^rn\dt + c{x'^ ,uj)dt, uj ^ fi. 

It is then natural to consider the following problem, where m is replaced by the proper 
measure P: on a filtered probability space {Q,J^,J^t,B,xo,Q), endowed with a Brownian 
motion B and with a J-q measurable random variable xq, we introduce the following system: 

ft T_ r __fjj n 1 _i _ 1 rt 



X7 



xo + /o 6K, Ps] ds + /o c(x^, w) ds + Bt, 



(22) ^ (J ~ /i, 

Pt = £(xt,w),vtG [o,r]. 

Proposition 4.2. 

There is pathwise existence and uniqueness for Equation ([22 



Proof. The proof is the same as given in Sznitman [26], Th 1.1, p. 172, up to minor modi- 
fications. The main idea consists in using a Picard iteration in the space of probabilities 
on C([0, r],S^ X R) endowed with an appropriate Wasserstein metric. We refer to it for 
details. El 

4.3. Fluctuations in the quenched model. The key argument of the proof is to explicit 
the speed of convergence as A^ — t- cxd for the rotators to the non-linear process (see Prop. 

KM- 

A major difference between this work and [12] is that, since in our quenched model, 
we only integrate w.r.t. oscillators and not w.r.t. the disorder, one has to deal with 
remaining terms, see Z^v in Proposition 14.31 to compare with [T2], Lemma 3.2, that would 
have disappeared in the averaged model. The main technical difficulty of Proposition 
14.31 is to control the asymptotic behaviour of such terms, see (f2^ . As in [12], having 
proved Prop. 14.31 the key argument of the proof is a uniform estimation of the norm 
of the process r/^'^^-*, see Propositions 14.41 and H^ based on the generalized stochastic 
differential equation verified by ry^'"^-*, see (|3Up . 

4.3.1. Preliminary results. We consider here a fixed realization of the disorder (cj) = 
{uji,uj2, ■ ■ ■)■ On a common filtered probability space (0, T, Tt, {B^)i>i, Q), endowed with 
a sequence of i.i.d. J^^-adapted Brownian motions (Bi) and with a sequence of i.i.d. Tq 
measurable random variables (i^*) with law A, we define as j;*'^ the solution of (|lip . and 
as x"^* the solution of (|22p . with the same Brownian motion B^ and with the same initial 
value ^*. 

The main technical proposition, from which every norm estimation of rj '^^> follows is 
the following: 
Proposition 4.3. 



(23) E 



sup 

t<T 



i,N 



<C/N + Zn{uji,...,un), 
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where the random variable {uj) i— t- Z]\}[uj) is such that: 

(24) lim limsupP (iVZAr(a;) > A) = 0. 

The (rather technical) proof of Proposition 14.31 is postponed to the end of the document 
(see f^ . Once again, we stress the fact that the term Zj\f would have disappeared in the 
averaged model. 

The first norm estimation of the process t] '^^' (which will be used to prove tightness) 
is a direct consequence of Proposition 14.31 and of a Hilbertian argument: 

Proposition 4.4. 



Under the hypothesis {H^ ] on /x, the process tj^'^'^' satisfies the following property: for all 
T>0, 



(25) 
where 



supE 

t<T 



Vt 



N,{u}) 



'3,2a 



< An{u:i,...,ujn), 



lim limsupP(^Ar > A) = Q. 



3,2o 



Proof. For all ip G Wg' , writing 



N,(u) 

Vt > V5 



1 ^ ■ 1 ^ 



--■■St [^)+Tt '{ip), 



we have: 

(26) 

But, by convexity, 



r]t ' \ipj <2[^St ' \ip) +T^ '(ip) j . 



N 



S^'^-\^f<Y^Dl. 



4 = 1 



Then, applying the latter equation to an orthonormal system ((^p)p>i in the Hilbert space 
Wq ' ", summing on p, we have by Parseval's identity on the continuous functional D ,:,jv ^i , 



E 



^N,{uj) 



-3,2a 



< E 



N 



(27) 
(28) 



"^t I'^t 1^'i- 



y\\D.. 

/ > II Xj 
i=\ 

N 

<C7^(l + |a;i|^")E 

i=l 

N 

< C J^ (1 + |wi|4") {C/N + Zn{uju . . . , oon)) , 



-3,2a 



i,N LOi 



i=X 



where we used (fTOll in (1271). and (123)1 in (1281). 
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On the other hand, 



E 



r,^'(-V)^ 



N 
1 



-E 



N 



E 



+ 


1e 


< 


> 


< 


> 



U=i 

AT 

,i=l 

^((/?(x^%a;i) - (Pt , <f)){ip{xt',ujj) - {Pt , if)) 



N 



Y,{^{xT^^i)^ + {Pt^vy 



1=1 

N 



+ l^Y.G{v)^^i)G{^){^,), 



i¥=j 



Y^ip{Xt\UJif 



i=l 



+ 2{Pt,vy + 



1 ^ \^ 



where G{if){u;) := J ip{y, uJi)P^' ( dy) — {Pt , if). If we apply the same Hilbertian argument 

as for S '^'^' , we see 

(29) 



E 



-,N,{iv) 



-3,2a 



2C^ 
< — E 

- N 



N 



E(i + 



CO, 



|4o^ 



,j=l 



+ C + 



1 ^ 

j=l 



-3,2a 



It is easy to see that the last term in ([^ can be reformulated as Bn{uji, . . . ,ujn), with 
the property that limA^oo hm sup^v^o^ P(-BAr > ^) = 0. Combining ([MD, dMD, ^B and 
(p9|) . Proposition 14.41 is proved. D 



4.3.2. Tightness of the fluctuations process. Applying Ito's formula to (jlip . we obtain, 
for all If bounded function on S^ x R, with two bounded derivatives w.r.t. x, for every 
sequence (w), for all t < T: 



(30) 



Af,(aj) \ / N,(lj) 

Vt , f) = iVo , f 



r/f'(-^/:fM)d. + Mf•(-)(^), 



where, for all y G S^, vr G R, 



C: {ip){y, n) = -^"{y, vr) + ^'{y, vr) {b[y, i/f ] + c{y, vr)) + (P, , (^'(•, ■)b{;y, vr)) , 
and Mj (c^) is a real continuous martingale with quadratic variation process 

M^'H(c^))^ = j[|'(.f'M,v.'(y,vr)2)d.. 

Lemma 4.5. 

For every N , the operator C^ defines a linear mapping from S into S and for all (p € S, 



^fi^) 



3,2a 



<c\\^\ 



6, a ■ 
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Proof. The terms ^ip"{y,TT) and ip'{y,iT)b[y,i/^] clearly satisfy the lemma. We study the 
two remaining terms: 



{Ps,^'b{.,y,7r))\\l = Y: 



ki+k2<3 



Ps, v'd k^d^k2b{-,y,Tr) 



S^xR 
1 



1 + k 



4a 



dydvr, 



^^ I 1^; |4a d^/ ^'iy,7rfPs{dy,d7T), 

J-R^ + m JsixR 

^^Il^llc3,. / ,J ua '^^ f (l + |vrr)2p,(dy,d7r), 



<C'II^I 



JR J^ + Fl JR 



And, 



9''(y>7r)c(y,7r) 



1 3,2a 



E 

fcl+fc2<3 



d k^d^k2 {ip'{y,TT)c{y,Tr)} 



S^xR 



1 + k 



4a 



■ dy dn. 



It suffices to estimate, for every differential operator Di = dy^t dj^^i , i = 1,2 with U1 + U2 + 
^1 + ^2 < 3, the following term: 



L 



SixR 



\Dnp'{y,7r)D2c{y,7r) 

1 + kl^" 



dy dvr < 



<C\\ip 



|Z)iy^(y,7r)nZ)2c(y,7r)|2(l + k '"^^ 
S^xR (1 + kl")' l + kl^° 

supyesi l^2c(y,7r)|2 



■ dydvr, 



1 6, a 



R 



1 + kP" 



dvr. 



The result follows from the assumptions made on c. 



D 



For the tightness criterion used below, we need to ensure that the trajectories of the 
fluctuations process are almost surely continuous: in that purpose, we need some more 
precise evaluations than in Prop. 14.41 
Proposition 4.6. 

The process (M^ ) satisfies, for every (uj), and for every T > 0, 



E 



sup 

t<T 



M 



N,{uj) 



-3,2a 



„ N 



to, 



l4a\ 



1=1 



Remark 4-^ ■ 

In particular, a consequence of (H^) is that, for P-almost every sequence (w). 



(31) 



supE 

N 



sup 

t<T 



M, 



N,{u;) 



-3,2a 



c 



N 
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Proof. Let {^p)p>i a complete orthonormal system in Wq' ". For fixed N, by Doob's 
inequality, X]p>iE \sup^^J'{M^ {^p))"^] is bounded by 



p>l 

N 



IVe 






(x:'^,w,)'ds 



A^ 



C 



i=l 

N 



r r^ 




2 




H ..iv 


ds 


Jo 


2-'s .^i 


3,2a 



^^E(i + i'^*i'")'(^si^g(izi 



N 



D 



i=l 



Proposition 4.8. 

For every N, every (uj), 

(32) 
with 



E 



sup 

t<T 


N,{u) 


2 

-6,Q 



< Cn{u>i, . . . ,u}n), 



lim limsupP(C7v > ^) = 0. 



Proof. Let (ipp) be a complete orthonormal system in Wq'°' of C°° function on S"*^ x R with 
compact support. We prove the stronger result: 



E 



^sup(77j , V'p 



p>i 



t<T 



< oo. 



Indeed, 



By Doob's inequality, 



E 



p>i 



t<T 



^sup/ry^ , V'p 



<C7 E 






P>1 



+ 5;e[m^'(-)(^,) 



p>i 
By Lemma 14.51 we have: 

Then, 

rT 



< C 



'Is 



-3,2a 



V'll6,a- 



E 



P>1 



<CM E 



<C7^T sup E 



Af,(w) 



Ar,(w) 



-3,2a 

2 
-3,2a 



ds. 



< C^TAn, 
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where A^r is defined in Proposition I4.4i The result follows. D 

Proposition 4.9. (1) For every N , forV-almost every (w), the trajectories of the fluc- 
tuations process rj^'^'^i are almost surely continuous in S' , 
(2) For every N, forF-almost every (lo), the trajectories of M^'^'^> are almost surely 
continuous in S' . 

Proof. We only prove for M^'^'^\ since, using Proposition 14.81 the proof is the same for 
^N,{ui) ^ Let {ipp) be a complete orthonormal system in Wq ' °, then for every fixed N 
and (a;), we know from the proof of Proposition 14.61 that for all e > 0, there exists some 
Mo > such that 

5^ sup(Mf(-)(^,))2 <£,«... 
Let {tm) be a sequence in [0, T] such that tm -^m-^oo t. 






p>i 



Mo 



p=i 

if tfn is sufficiently large. D 



We are now in position to prove the tightness of the fluctuations process. Let us recall 
some notations: for fixed A^ and (w) Ji^^y^) is the law of the process ry-^'V'^). Hence, Ji^^y^) 
is an element of M.i{C{[0,T],S')), endowed with the topology of weak convergence and 
with B*, the smallest a-algebra such that the evaluations Q ^ {Q , f) are measurable, / 
being measurable and bounded. 

We will denote by Qn the law of the random variable (w) i— )• 'H^'^'^' . The main result 
of this part is the following: 

Theorem 4.10. (1) forW-almost every sequence (cj), the law of the process M^'^-^' is 
tight inMi{C{[0,T],S')), 
(2) The law of the sequence (w) ^ ?^^'M is tight on Mi {Mi{C{[Q,T],S'))). 



Before proving Theorem 14. 10^ we recall the following result and notations (cf. Mitoma 
[19], Th 3.1, p. 993): 
Proposition 4.11 (Mitoma's criterion). 

Let {Pn) be a sequence of probability measures on (Cs' := C{[0,T],S'),Bcg,)- For each 
ip & S, we denote by H^ the mapping of Cg' to C := C([0,T],R) defined by 

n^ : V'(-) e Cs> ^ (V(-) ,v)ec. 

Then, if for all ip £ S, the sequence (PatH"^) is tight in C, the sequence {Pn) is tight in 
Cs'. 

Remark ^.12. 

A closer look to the proof of Mitoma shows that it suffices to verify the tightness of 
(P/vH"^) for c/3 in a countable dense subset of the nuclear Frechet space (5, || • || ,p > 1). 
Thanks to Mitoma's result, it suffices to have a tightness criterion in R. We recall here 
the usual result (cf. Billingsley 0]): A sequence of (0^,7"/^)-adapted processes {Y^) with 
paths in C([0, T], R) is tight if both of the following conditions hold: 

• Condition [T]: for all t < T and (5 > 0, there exists C > such that 

^Ac) supP(|y,^|>C7)<5, 

N 
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Condition [A]: for all r]i,rj2 > 0, there exists C > and Nq such that for all 
J^^-stopping times tn, 



i^vi, m,c) 



sup SUpP(|y,^-y,^+e|>7?2)<7?i. 
N>No 9<C 



Proof of Theorem \4.10[ (1) Tightness of {M '^^'): for a fixed realization of the disor- 
der (w), for fixed ip £ S, we have: 

• For ah t £ [0,T], for all 6 > 0, for all C > 0, 



Mf'(-)(^) 



>C] < 



E 



SUPi<T 



{Mf(-)M2} 



C2 



E 



< 



SUPi<T 



m; 



Af,(a;) 



I-3.2Q 



l|2 
'/'113.2a 



C2 



/^ II IM 
^ f-- IIV5|l3.2a 

- 2r~~ ^^^ AT 



1 ^ 

-j;(i+kr),(cf. m), 



i=l 



for a suitable C > (depending on (f^i^)). Condition [T] is proved. 
• Let us verify Condition [A]: For every ip £ S, for every 6,9,rii,ri2 > 0,0 < 6, 
for every stopping time tn, 



m 






TJV- 



{^s , v'{y,T^f) ds 



TN 

1 



< (O ^ — E 
— Wr' 116, a „2 

'/2 L-J TN 



TN- 



'/2 



n ,,2 (^(J 
^ IIV'll6,a-;2-SUP 
V2 N 



S^xR 



Af 



^1^,77 II %„ di/f ds 






™ " i.i 



Wi 



|4o\ 



ds 



, (cf. (HBD and (HH)), 



(^E(i + N*")) 



This last term is lower or equal than r]i for 6 sufficiently small (depending on 

(OJ)). 

(2) Tightness of (Qn)'- we need to be more careful here, since the tightness is in law 
w.r.t. the disorder. Let {^j)j>i be a countable family in the nuclear Frechet space 
5. Without any restriction, we can always suppose that {{(pj ||g^ = 1, for every 
j > 1. We define the following decreasing sequences (indexed by J > 1) of subsets 
ofA^i(C([0,r],5')): 

Kliipi, ...,^j):=[P; Vt, VI <j<J, PU-j satisfies (Tt,5,c^)} , 
K^{ipi,...,(pj) := |P; VI < J < J,Vr/i,??2 > 0,PU~j satisfies (A,,i,,,2,C2)| , 

where Ci = Ci(e,(5), C2 = C2(e, 771,7/2) will be precised later. By construction and 
by Mitoma's theorem (cf. Remark 14. 12p . 
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is a relatively compact subset of TMi (C([0, T],5')). In order to prove tightness of 
{@n), it is sufficient to prove that, for all e > 0, 



(33) 



Vi = 1,2, limsupeAT \[\Kf{ 
^ VY 



h,...,ct>jr\ <e. 

For e > 0,let A = A{e) such that lim inf Ar_^oo P {An < ^) > 1 - e, and 
liminf P ( 1 V (1 + \uJi\^n + ANicoi, . . . , wjv) < ^ ) > 1 - e, 

\ i=l / 

where A]\f is the random variable defined in Proposition 14. 4i We define the corre- 
sponding constants (for a sufficiently large constant C): 



Ciis,5):-- 



A{e) 



C2i£,m,m) 



Then, 



/ 



QNiKli<Pi,...,<Pj))> 



E 



(u), vt, VI < i < J, ys, 



m^2 

CA{e) ' 



n^'^' , 0. 



V 



Cii5,e) 



<d 



J 



>P (uj), supE 



t<T 



Vt 



N,{uj) 



-6,a 



< A\ , (by definition of Ci ) , 



>F{An < A{e)) , (cf. (HE]) and ([25])). 

Letting J — ;■ oo in the latter inequality, we obtain: Qj\f{\Jj Kf^cpi, . . . ^cpjY) < 
F{A]\f > A). Taking on both sides limsupjv^o^, we get the result. 

Furthermore, for r]2 > 0, < < C2 and t^ < T a stopping time, for all 
1<J<J, 



TN+d 



r,^'(-\ cf i^,)) ds 



\ 


1 


i-TN+e 


V2 
J 


'r 


Jtn 






Jtn 



r,^'(-\ cf i^,)) ds 



v^'^-\cf{^,) 



ds 



m Jo 

CCo '^ 



v^'^-Kcfi^,) 



ds. 



< 



V2 Jo 



E 



N ||2 



ds, 



< 



CTC2 



An, (cf. (USD). 



And, 



P (|<+.(^.) - M?^{^,)\ > ^2) < ^ (^^ E (1 + l-^l'")) • 



So, for all J > 1, by definition of C2, 



^fA1(v^.)-.f.(-)(^.) 



>r/2l < 



A{e] 



(^^iv + ^E(i + i'^^i'"))- 
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Consequently, 



( 1 ^ 

Qf,{Kl{^x, ..., fj)) > IP Utv + ^ j; (1 + 

V i=i 



UJ; 



|4o\ 



Letting J ^- cc, we get limsup^Y Qn (Uj -^llv'i) ■ • • ' '/'j) 



< e. Eq. ([33|) is proved. 

D 



4.3.3. Identification of the limit. The proof of the fluctuations result will be complete when 

we identify any possible limit. 

Proposition 4.13 (Identification of the initial value). 

The random variable (uj) i— )• £ ( t^q I converges in law to the random variable uj i— )> 

C{X{lo)), where for all oj, X(uj) = C{uj) + Y , with Y a centered Gaussian process with 
covariance Fi. Moreover uj i— t- C{u}) is a Gaussian process with covariance T2, where Fi 
and F2 are defined in ([8|) and ([9]). 



N,{uj) 



N,(uj) 



N.(u)) 



Proof. For simplicity, we only identify here the law of ( ?7g , ^) for all if. The same proof 

works for the law of finite-dimensional distributions 
p > 1. We write Fj for Fj(c^, ip), i = 1,2. One has: 

N 

ip{x,u!i)X{dx) 



N,(ui) 

Vo ,f 



1 ^ 



i=l 



ip{x,u}i)X{dx) - {uo, If) 



It is easy to see that B^'^'^' converges in law to Z2 ~ AA(0, F2). Moreover, for P-almost 

every (w), A^'^"^^ converges in law to Zi ~ A/'(0,Fi) (see Billingsley [5], Th. 27.3 p. 362). 

2 

That means that for all u G R, '0a^(^) •= ^A (e*" ' j converges to V'y(^) '■= ^ ^""^ • 
But, then, for all F € Cfe(R), 



E 



F{ E 






E 



F E 



'-'X 



„iM(A^'(")+B'V'('^)) 



E 



F(e»^"''"V^(n) 



Since ^n{u) converges almost surely to a constant, the limit of the expression above exists 



(Slutsky's theorem) and is equal to E 



F[e 



iuZ2-^ 



n 



Proposition 4.14 (Identification of tfie martingale part). 

For ¥ -almost every [uj], the sequence (M^'''^-*) converges in law inC{[0,T],S') to a Gauss- 
ian process W with covariance defined in ([7]) . 

Proof. For fixed (uj), [M '^^') is a sequence of uniformly square-integrable continuous 
martingales (cf. Remark 14. 7|) . which is tight in C([0,T],5'). Let Wi and W2 be two 
accumulation points (continuous square-integrable martingales which a priori depend on 
(cj)) and (M'^(^)'('^)) and (M'^(^)'('^)) be two subsequences converging to Wi and W2, 
respectively. Note that we can suppose that </)(A^) < ip{N) for all A^. For all ip £ S, 
lim^^oo (M'^W'H(v?) , M^(^)'('^)((^))^ = {Wi{(p) , W2{(p))t, for all t, and 



/^0(iV),H(^)^^^(7V),H(^) 



A^) 



Av' 



/\2 



ds. 
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We now have to identify the hmit: we already know that for P- almost every realization 
of the disorder (w), (v '^'^>) converges in law to P. But, the latter expression, seen as 
a function of u, is continuous. So {Wi{ip) , W2{^))i = /q {Ps , i^')^)- So Wi — W2 is a 
continuous square integrable martingale whose Doob- Meyer process is 0. So Wi = W2 and 
is characterized as the Gaussian process with covariance given in ([7]). The convergence 
follows. D 

Proof of the independence of W and X. We prove more : the triple (Y, C, W) is indepen- 
dent. For sake of simplicity, we only consider the case of {Y{^p),C{ip),Wt{^)) for fixed t 
and If. 

Let us first recall some notations: let A^'^'^', B^'^'^' and M^ (97) be the random 
variables defined in the proof of Proposition l4.13l and l4.14l and let ■0a^ (^) •= E ( e*""^ " ) > 

i>B^{^) ■= E fe™-^'^' '^' j, VAf^(^) '■= E (e™^^t '" ("^M be their characteristic functions 

{u,v,w G R). We know that, for almost every (w), iP^n{u) converges to ipriu) = e ^'"i 
and that ipj^^N^w) converges to the deterministic function il^wiw) := E (e™^*('^-'). But, if 
i^civ) = IE (e*""-^), then, for all u,v,w E R, using the independence of the Brownian with 
the initial conditions. 



Using Slutsky's theorem, we see that any limit couple (Y, C, W) satisfies 

E (e (^e'^'^+i^c+iwWti^)'^'^ ^ 4;Y{u)iPc{v)iJw{w). 
which is the desired result. D 

We recall that the limit second order differential operator Cs is defined by 

^s{f){y,T^) ■■= -f"iy,TT)+ip'iy,7r){b[y,Ps]+c{y,TT)) + {Ps , ip' {■ , ■)b{- , y , ir)) . 

As in Lemma 14.51 we can prove the following: 
Lemma 4.15. 



,N 



Assume (Hb^c)- Then for every N , s < T, (co), the operator Cs and C^ are linear 
continuous from S to S and 

II Ce(ijy) lie < C II CO llo 

II '-s\Y'J 116, a — ^ II r 118, a ' 



,N , 



We are now in position to prove Theorem 12.101 

Proof of Theorem \2.1(K Let G be an accumulation point of Qn- Thus, for a certain sub- 
sequence (which will be also denoted as N for notations purpose), the random variable 
(w) I—)- Ji"''^'^) converges in law to a random variable H with values in A^i(C([0, T],S')) with 
law ©. Applying Skorohod's representation theorem, there exists some probability space 
(Q^^), P^-*^', J^'^)) and random variables defined on Q^^', ioi 1— ?■ H^{iOi) and coi 1— t- H{loi) 
such that H^ has the same law as (w) 1— )• 'H^'^'^', H has the same law as Ti, and for 
p(i)-almost every ui G f](i), H^{ui) converges to H{uji) in Mi{C{[0,T],S')). 

An easy application of Proposition 14.81 and Borel-Cantelli's Lemma shows that P'^)- 
almost surely, E ( supj<y || 77^^ ll-6a) ^ °*^- Then we know from Lemma 14.151 that the 

integral term L C*r]^^ ds makes sense as a Bochner's integral in Wq '" C S' . 

Let r/^'^i with law H^ {uji); rj^''^^ converges in law to some rj'^^ with law HiuJi). By 
uniqueness in law convergence, using Propositions 14.131 and I4.14| we see that (t^q ^ , W) as 
the same law as {X{uji),W). For fixed (/? G 5, we define F^p from C([0,r],cS') into R by 
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Fipil) '■= {it ) ^) — (70 ^ v) ~ Jq {is ) J-'sV') ds. The function F^^ is continuous and since 
^N,u}i converges in law to r/'^^, the sequence (-^(^(?? )) converges in law to F^(r/'^^). To 

prove the theorem, it remains to show that the law of the term L ( rjs '^^ , C^ 93 — Cs^p ) ds 
converges in law to 0. We show that there is convergence in probability: For all e > 0, for 
all A> 0, using Proposition I4.8| Lemma [4. 151 and Cauchy-Schwarz's inequality, 

Un,s ■■= P^'^ f E ' ' 
E 






ds 



> e 



> e 



< 



< 



< 



E 



JV,(w) 



-6,Q' 



1/2 



E 



CN{i^l,...,iOMf'^ j E 



E 



(^^ -^s)M 



1/2 

ds > 



6, a 



1/2 



ds > e 



2 

6, a 



1/2 



A 



ds > e (cf. Prop SS 

+ P(C^>A). 



Using (14. Sp . it suffices to prove that, for all e > 0, 



(34) 



lim sup ] 

7V-5>oo 



E 



(^^ -^.)M 



6, a 



1/2 



ds > £ 



0. 



Indeed, for every ip £ S, 

UfiifMn) := (cf -Cs){^){y,7r)=^'{y,7r){b[y,u^]-b[y,P,]). 

An analogous calculation as in Lemma 14.51 shows that, using Lipschitz assumptions on b, 
and Proposition 14.31 

2 



E 



snpU^iip) 

s<t 



6,0 



< \\y^\\'{C/N + Dn{uji, . . . ,un)), 



with the property that lim^_j.oo lim sup^ F{NDn > ^) = 0. Equation (fH^ is a direct 
consequence. 

Since there is uniqueness in law in ([TU]) . is perfectly defined, and thus, unique. The 
convergence follows. D 

5. Proofs for the fluctuations of the order parameters 
We end by the proofs of paragraph 12.3.31 

5.1. Proof of Proposition 12.131 

(1) This is straightforward since r '^^> = ^l^u '^^' , e*^)| and since for P-almost every 
disorder (w), i/^'('^) converges weakly to P. 

(2) The following sequences are well defined: VA; > 0, 



Uk{t):= f e-l^lw^cos(0)Pt(d0, dw), 

JSixR 
Js^xR 



LetS=(^oo(N), 
with its usual 11 • 11 



iioo 

norm 



be the Banach space of real bounded sequences endowed 



supfc>o l^fcD- For all t > 0, let .At : ii^ x £' 
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Ex E,he the following linear operator (where (u, v) is a typical element oi Ex E): 
For all /c > 

^(-",0)^ = -^Uk - ak{t)vQ + Pk{t)uQ - Kvk+i, 
A(0,v)fc = -i2Vk + '~lk{t)vQ- ak{t)uQ + Kuk+i, 



where, 



Qfc(t) = (Pi,e-H^fccos(-)sin(-)) 
/3fc(t) = (Pi,e-H^^sin2(.)), 
7fc(t) = (Pt,e-H^'=cos2(.)). 



(t, u, v) I— >• ^i • (u, v) is globally Lipschitz-continuous map from [0, T] x E x E into 
E X E and one easily verifies considering ([!]) (in the case of the sine-model) and 
developing the sine interaction that 1 1— )■ (u{t),v{t)) satisfies in E x E the following 
linear non-homogeneous Cauchy Problem: 



Uk{0) 
Vk{0) 



Af{uit),v{t)) 



w'^cosl 



Pn 



-\ui\l ,k 



uj sm • 



V/c> 
\/k> 0. 



Let us suppose that there exists some to £ [0, T] such that rj^ = 0, namely uo(to) = 
^o(^o) = 0. Then, if (n, v) is the constant function on [0, T] such that for all /c > 0, 
% = Uk{to),Vk = Vk(to), then {u,v) satisfy the same Cauchy Problem as {u,v) 
with initial condition at time to- By Cauchy-Lipschitz theorem, both functions 
coincide on [0,T]. In particular, uq and vq are always zero and thus r = 0. 
(3) We suppose ^H^. A simple calculation shows that the fluctuations process TZ^ 
verifies for all t £ [0, T], 



n 



N,{ui) 



Vt 



N,{ui) 



, cos(-) 



NXoj) 



+ Pi,cos(•)) + (r?f•^"^sin(•) 



N,{ui) 



+ Pt , sin(-) 



N,{Ld) 



+ rt 



3? 



^t 



N,{u)) 



N,(ui) 



+ Pt 



N,{co) 



, e' 



+ n 



Let u^'('^) := (i/^'(^), e*^), u^'(^) := (r/^'('^) , e^^) and u := (P,e^^), u"' : = 
^^o; ^ g«x^ ijg their corresponding limits. The result follows if we prove the follow- 



ing property: the random variables (w) i— )• £ (n 



N,(u,)^^N,(u,) 



)) converges in law to the 



random variable uj >-^ C{u, v'^). The tightness of this random variable follows from 
the convergence of both empirical measure and fluctuations process. As already 
said in Remark 12.141 it suffices to prove the convergence of the finite-dimensional 



N,{ui) N,{uj)^ 



N,{^) 



,---,Ut 



N,{uj)s , N,{ui) 






N,{u;). 



) • • • ) '^i^ 



for all ele- 



marginals (n^ ,Vf- ) = f(^tti' 
ment of [0, T], ti, . . . ,tp, p > 1. 

Since the limit of (u^''^^) is a constant, this is mainly a consequence of Slutsky's 
theorem. But since this is a convergence in law with respect to the disorder, one has 
to adapt the proof. We prove the following: VG G C^(R), Vr = (ri, . . . ,rp) G R^, 
Vs = (si,... ,Sp) £ W, 



E 



G( V'(^iv,(-)_^iv,(^)-j(r,s) 



E 



N-^oo 



G(93(„^,^^)(r,s)j 
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where ^[x,Y)(—^^ ~ ^ [e*-— "*"*-— ] is the characteristic function of the couple 
(X, y). Indeed, we have successively: 



a-N ■-- 



< 



+ 



E 
E 
E 



^('/'(„f("),.f-("))fc^) 



G{ ff^^N,(u,)^^N,(u.)^{r,s) 



-E 
E 



G[^iu„vf){r,s) 



^('^fe,.r'(-')(^'^) 



E 






G[V(u„vi'){r,s) 



< CE 
E 



'^(nf-("',.f'(-))(^'^) - '^{n„.^("')fc^) 



+ 

< CEE 



G ip N.{c.)Ar,s) 



E 



N,(uj) 

e — i- — e — i- 



E 



^('^fe,.f("')fc^) 



E 



G{v{u,_,tl){L:S)) 



G[^{u,_,v^){L:i 



But, we have E 



e -i -i — e — ^^ 



< min ( 2, |r||u^ — Mtl ) • So, for all e > 0, 



E 



NAuj) 

e — i- — e — ^ 



<e\r\ +2P 



Ar,(w) 



Uj 



> e 



Taking lim sup^_j.oo , and letting e — )• 0, we get lim a^r = 0. The result follows. 



5.2. Proof of Proposition 12.151 The proof is similar to the previous one and relies on 
the two following equalities : 



c 



N,(u]) 






^ (-;"'"' - () = 777^ (^ (""''"' ■ '") + <^ ■ '"> '^"'"0 ■ 



Appendix A. Proof of Proposition 14.31 

Thanks to the Lipschitz continuity of b and c, introducing u as the empirical measure 
corresponding to {x'^\uji), we have, (inserting 6[a;^% i^^] — b[x'^',Us] in the b term). 



E 



sup|xi'^-x^f 

.s<t 



+ / E 

Jo 



b[xi'\i^^]-b[xT,Ps]y 



ds 



c(x:'^,^,)-c(x^% ..,))'■ 



ds], 



+ / E 

JO 



sup|x^^-x-f 



{b[x'^^,u,]-b[xT,Ps]r 



ds + / sup E 

'o l<j<N 

ds] . 



SUp|x„ -X:^' I 
u<s 



ds 
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Applying Gronwall's Lemma to supx<j< jv E 
for some Zjy: 



sup„<i 



t^i j,N 

u ■^u 



it suffices to prove that 



Jo 



E 



{b[x'^\us] - b[x^\Ps]f ds < C/N + Zn{uoi,...,uon)- 



Indeed, for all 1 < i < A^, (we write x* instead of x^^ to simplify notations): 



,2 1 



N 



u,,N ■•= {b[xl,Us]-blxi,Ps^' = -2 Y.^{x\x^f + Y,nx\x'')T{x\: 



vJ=l 



k^l 



where T{x^,x^) := b{x\,xi,ujj) — Jb{xl,y,7r)Ps{dy, dvr). Since b is bounded, we see that 
the first term is of order (1/iV). We only have to study the remaining term: 



E 



^T(x^x'=)^(x^x') 



kf^l 



<CN + E 



Y, T{x\x^)T{x\x^) 






Since the (x*) are independent, if we take conditional expectation w.r.t. {x^',r ^ I) in the 
last term, we get: 



E 



T{x\x'')T{x\x^) 






E[E[Y^T{x\x'')Tix\x': 



E 



x^r// 







" 


)_J T{x\x'^)Gi{x') 


= E 


> ; Gk{x^)Gdx^) 


k^i,l^i 




kj^i.ly^i 



where Gi{x) = G{x,uji) = fb{x,y,ui)P^'{dy) - fb{x,y,TT)Ps{dy, dvr). Defining 



ZNiiOu...,u^N):=^£B (-l=f2G{xl,ui)\ 



ds, 



in order to prove ()24p it suffices to show that for some constant C, 



E 



E 



1 ^ 

1=1 / 



ds 



< C. 
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The rest of the proof is devoted to prove this last assertion: we have successively (setting 



E 



Jo 



Nyx'^',^^) ds 



Jo 

' N N 

fc=l 1=1 

' N 



^l> 


E 
J 


^I> 


E 



ds, 



1=1 



ds + C 



H / E 



E 



l^k, l^i, k^i 



ds. 



The first term of the RHS of the last inequality is bounded, since b is bounded. But, if 
we condition w.r.t. uJr for r ^ i,r ^ k, we see that the second term is zero. The result 
follows. 
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